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ABSTRACT 

The mass at which a transition is made between stars that have radiative or 
convective cores throughout the core H-burning phase is a fairly sensitive function 
of Z (particularly, the CNO abundances). As a consequence, the ~ 4 Gyr, open 
cluster M 67 provides a constraint on Z & (and the solar heavy-element mixture) 
because (i) high-resolution spectroscopy indicates that this system has virtually 
the same metal abundances as the Sun, and (ii) its turnoff stars have masses 
just above the lower limit for sustained core convection on the main sequence. 
In this study, evolutionary tracks and isochrones using the latest MARCS model 
atmospheres as boundary conditions have been computed for O.6-I.4.M0 on the 
assumption of a metals mix (implying Z & pa 0.0125) based on the solar abun- 
dances derived by M. Asplund and collaborators using 3-D model atmospheres. 
These calculations do not predict a turnoff gap where one is observed in M67. 
No such difficulty is found if the analysis uses isochrones for Z Q = 0.0165, assum- 
ing the Grevesse & Sauval (1998) mix of heavy elements. Our findings, like the 
inferences from helioseismology, indicate a problem with the Asplund et al. abun- 
dances. However, it is possible that low-Z models with diffusive processes taken 
into account will be less problematic. 

Subject headings: Hertzsprung- Russell diagram — open clusters: individual (M 67) 
— stars: atmospheres — stars: evolution — Sun: abundances 
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Introduction 



The metallicity of the Sun (both the mix of heavy elements, and the total mass-fraction 
abundance, Z) is presently a subject of considerable controversy. Until a few years ago, the 
generally accepted value of Z Q was ~ 0.017-0.02, and there was good consistency between 
the predictions of Standard Solar Models (SSMs) for meta llicities within this range (see, e.g., 
Turcotte et al.lll998l . iBahcall. Pinsonneault. fc Basull200ll ) and the measured neutrino flux as 
well as helioseismological data. However, significantly reduced abundances for several of the 
most abundant heavy elements (including C, N, O, and Ne), resulting in Z & ~ 0.0125, have 
been derived by M. Asplund and collaborators from their analyses of the photospheric spec- 
trum using 3-D, non-LTE model atmospheres (see I Asplund. Grevesse. fc Sauvall 120061 ; and 
references therein). Whereas high-Z solar models are able to reproduce the inferred radial 
dependence of the square of the sound s peed (down to ~ 0.1R Q ) from solar oscillations to 
^ 0.3% (see lChristensen-Dalsgaardll2002l ). thos e computed for the Asplun d et al. metallicity 
are unable to do so to better than ~ 3% (e.g., (ITurck-Chieze et al.l 12004 ). 



It seems unlikely that this difficulty is due to problems with current opaci ties (lAntia fc Basu 



2005 f) or the presently accepted rates of diffusive processes (IGuzik! 120051 ) . As noted by 
Montalban et al.l (120061 ). any substantial increases to the diffusion coefficients that are made 
in order to recover good agreement with the inferred sound speed profile result in predictions 
for the surface helium content of the Sun that are below the values obtained by the inversion 
of solar oscillations. In fact, turbulent mixing processes appear to operate in stellar envelopes 
to reduce the effects of diffusion on the surface abundances of ~ lAi & main-sequence (MS) 



stars (see 



Korn et al.ll2006l : iRichard et al J 12002( 1 . 



While the low solar Z derived by Asplund et al. may be erroneous, the 3-D model at- 
mospheres used in their stu dies have had u nprecedented success in modelling the properties 
of the solar atmosphere (see I Asplund! 120051 ) . and the predicted spec tral line profiles base d on 
these atmospheres provide superb matches to the observed profiles (lAsplund et al.ll2000l ). In 
addition, the reduced abundances of the C NO elements are in excelle nt agreement with the 
values measured i n the interstellar med ium (ITurck-Chieze et al.ll2004l ) and in solar neighbor- 
hood B stars (see lAsplund et al.ll2005t and references therein). Thus, their findings are not 
easily dismissed. 

There is another, albeit indirect, way of constraining the metallicity of the Sun that has 
not yet been investigated. The large reduction in the CNO and Ne abundances found by 
Asplund et al. (relative to previous estimates) affects the mass at which a transition is made 
between stars that have radiative cores during the MS phase and those possessing convective 
cores when their central hydrogen fuel is exhausted. Consequently, a cluster like M 67, which 
has [m/H] « 0.0 and turnoff stars with masses close to the transition mass, A4 tT , provides 
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a probe of both stellar physic s and the chemical composition of the material out of whic h 
the stars formed. As shown by lVandenBerg. Bergbusch. &: Dowlerl (12006; hereafter VBD06I ). 
M. tr is predicted to increase with decreasing metal abundance, and to be a fairly sensitive 
function of Z (see their Fig. 2). Whether or not an old open cluster is expected to have 
a gap in its distribution of near-turnoff stars on the color-magnitude diagram (CMD) will 
therefore depend on its metallicity. (A gap is the observational manifestation of the rapid 
contraction phase that accompanies the exhaustion of hydrogen in stars that have convective 
cores at the end of their MS lifetimes.) 



The M 67 CMD (see Montgomery. Marschall. fc Janesl Il993l : ISandquistl |200J) contains 



such a gap, and stellar models for Z <^ 0.017 have been quite successful in reproducing 
these data. IVandenBerg Stetso n (200 41) hav e shown that a 4.0 Gyr isochrone for Z = 
0.0173, assuming the iGrevesse fc Noels! (119931 ) metal abundances, provides a satisfactory 
match to the cluster CMD, if a small amount of con vective core ove r shoot ing is assumed 
(also see VBD06). A nearly identical fit was obtained by lMichaud et al.l (120041 ) using diffusive 
isochrones for 3.7 Gyr and initial = 0.02, without having to postulate any core overshooting. 
These calculations predict that the cluster MS stars have ~ 6-12% surface underabundances 



of the metals, de pending on their tem peratures — compared with a 



< 



reduction in the 



case of the Sun (ITurcotte et al.lll998l ). which is much too small to simulate the difference 
between the Asplund et al. solar composition and that reported by, e.g., Grevesse & Noels. 

Given the apparent success of relatively metal-rich models in explaining the CMD of 
M67, which has very close to the solar metallicity, is it possible for Z = 0.0125 isochrones to 
provide a comparable fit? This is the question that has motivated the present investigation. 



Stellar Models for Solar Abundances 



The evolutionary tracks needed for this investigation have been computed using the 
Victoria stellar evolution code (see VBD06, and referen ces therein). Model a tmospheres 



produced by the latest version of the MARCS code (e.g., iGustafsson et al.l 120031 ) have been 
used as boundary conditions. [Differences in the treatment of the surface layers do not alter 
the tracks for the MS and early subgiant branch (SGB) phases of solar metallicity stars in 
any significant way if the value of the mixing- length paramet er, amlt, is chosen so th at the 
solar constraint is satisfied (see, e.g., Fig. 3 in the study by lVandenBerg et al.l 120071 ). The 
assumed value of «mlt does affect the predicted location of the red-giant branch (RGB), but 
that is not a concern for the present work.] 



The standard abundance distribution reported by 



Grevesse &: Sauval 



(11998: hereafter GS981 ). 
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as determined from analyses of solar photospheric spectra and meteoritic data using classical 
1-D hydrostatic models, was the last such compilation prior to the use of 3-D hydrodynam- 
ical model atmospheres by Asplund et al. (2006) to determine the solar composition. The 
GS98 abundances for the 19 metals normally considered in stellar evolutionary computations 
(because they are the only heavy elements that ar e taken into account in the calculation of 



OPAL opacities for stellar interior conditions; see llglesias fc Rogers! Il996l ) are listed in the 
second column of Table HJ The third column contains the the abundance distribution that is 
obtained when the revised values derived by Asplund et al. for several of the heavy elements 
are adopted: note that the C and N abundances in Table [1] are 0.02 dex higher than their 
published 2006 estimates. 

OPAL opacities were computed for both heavy-element mixtures using the Livermore 
Laboratory website facility (see |http: / / www-phys.llnl.gov/Research/OPAL ) . Because these 



calculations do not include the contributions from molecular sources or grains, complemen- 
tary opacity data for T < 1 4 K w ere computed for the same abundances using the code 



described by Ferguson et al.l (120051 ). The tabulated values of log Nn e were obtained from 
SSMs that were constructed for the two sets of abundances: the resultant values of Y and Z 
are listed in the last two rows of Table [TJ To satisfy the solar constraint, the 1.0A4 Q models 
for the GS98 and Asplund et al. abundances required Q!mlt = 1-84 and 1.80, respectively 



(whe n MARCS model atmospheres are used as boundary conditions; see IVandenBerg et al. 



20071). 



The value of Ai tr is a function of Z (especially the CNO abundances; see below). 
The transition occurs when the CNO-cycle becomes an important source of nuclear energy 
production, and since both the decrease in the abundances of the CNO elements and the 
concomitant reduction in opacity will serve to reduce the rate of the CNO-cycle, the minimum 
mass for sustained core convection on the MS is higher in stars of lower Z. By trial and 
error, A4 t r was found to be 1.195A^ Q if Z = 0.0125 (Asplund et al. metals mix), whereas 
1.155.M© is obtained if Z — 0.0165 (GS98 metal abundances). [The transition mass also 
depends on the adopted helium content; e.g., A4 tr = 1.181A4© if Z — 0.0125 and the value 
of Y required by an SSM for Z = 0.0165 is assumed (i.e., Y = 0.2676; see Tabled]).] 

Indeed, the character of the tracks changes abruptly when core convection persists until 
H exhaustion. This is shown in Figured! which plots the evolutionary tracks for 1.0— 1.4M© 
that were computed for the two values of Z Q ; including, in particular, sequences for A4 tr 
and (A4 tr + 0.001)A4 Q . Each of the higher mass tracks possesses a blueward hook, which 
arises when a star contracts following the depletion of hydrogen in a convective core in order 
to ignite H-burning in a shell around the He core. This morphology is not seen in the lower 
mass tracks because the central regions are radiative when hydrogen is exhausted and there 
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is a smooth transition to H shell burning. 

A difference of O.O4A^ in A4 tT is surprisingly large given that VBD06 found an increase 
of only 0.009A4 Q when Z was decreased from 0.0173 to 0.0125. However, the latter calcula- 



tions assumed the same relative abundances of the metals (those given by iGrevesse Sz Noels 



19931 ). whereas the tracks plotted in Fig. [T]assume two very different heavy-element mixtures. 
Indeed, a large reduction in the CNO abundances is mainly responsible for the downward 
revision of the solar value of Z from 0.0165 (GS98) to 0.0125 (Asplund et al.) - which 
demonstrates that A4 tr is a strong function of the abundances of the CNO elements. 

An increase in A4. tr from 1.155 to 1.195A4© will be accompanied by a decrease in the 
maximum age at which a gap near the turnoff is expected in an observed CMD. Fig. [1] shows 
that the biggest differences between the two sets of models occur at masses of ~ 1.2A4 Q . 
Because A4 tT is lower in the case of the GS98 computations, and because the amount of 
overshooting in stars of mass A4* > A4 tr is assumed to vary directly with Ai* — Aitr (see 
VBD06), the tracks represented by the dashed curves for masses of 1.2-1.4A4 have larger 
convective cores during MS evolution, and cooler blueward hooks at H exhaustion, than 
those plotted as solid curves. Such differences will impact the isochrones (computed in this 
study using the interpolation software described by VBD06) and the degree to which they 
are able to reproduce the M67 CMD. 



3. Application to the M 67 CMD 



Randich et al. 


2006L). For the cluster reddening. 


as found from the 


Schleeel. Finkbeiner. & Davis! 



0.0±0.03 fsee Tautvaisiene et al. 



0.038 



is 



mag, 

probably accurate to within ±0. 01 mag, given that it agrees so wel l with the values derive d 
using alternative methods (see iNissen. Twarog. fc Crawford! 119871 ; ISarajedini et al.lll999l ). 
As far as the distance is concerned, the modulus that is obtained from a main-sequence fit of 
the dereddened photometry to the lower MS segments of our isochrones should also be quite 
accurate (probably to within ±0.1 mag) because th e model temperatures satis fy the solar 
constraint and the adopted color transformations (by IVandenBerg fc Cleml 120031 ) are i n very 
good agreement with the empirical relations derived by ISekiguchi fc Fukugital (120001 ) from 
their study of solar neighborhood stars having well-determined properties (see the Vanden- 
Berg & Clem paper). The application of the MS-fitting technique yields (m — M)y = 9.70. 

With these choices for the basic parameters of M 67, the isochrone for each value of Z Q 
that provides the best fit to the subgiants is readily identified. As shown in Figure [2J the 
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3.9 Gyr isochrone for Z = 0.0165 (GS98 metal abundances) provides a good match to the 
morphology of the M 67 CMD all the way from the ZAMS to the base of the giant branch, 
including, in particular, the location of the gap near the turnoff. (Many of the stars just 
brighter than the gap, at My « 3.1, are presumably binaries, given that M67 is known 
to have an unusually large binarity fraction; at least 63%, according to Montgomery et 
al. 1993. The ~ 0.06 mag offset in color at a given My along the lower RGB is discussed by 
VandenBerg et al. 2007.) In stark contrast with this, the 4.2 Gyr isochrone for Z = 0.0125 
(the Asplund et al. metallicity) fails to reproduce the location of the near-turnoff stars of 
M67, despite providing a comparably good fit to the ZAMS and the SGB. This isochrone 
shows no indication at all of a blueward hook feature, and thus it does not predict a gap in 
the distribution of stars near the turnoff where one is observed^] 

The oldest isochrone in the Z = 0.0125 set that has the required turnoff morphology is 
one for an age of 3.6 Gyr, and indeed, it provides a satisfactory fit to the CMD of M 67 if 
E(B — V) ~ 0.065 and (m — M)y ~ 9.90. However, such a high reddening seems to be ruled 
out and it is doubtful that the model colors are in error by as much as 0.03 mag. Isochrones 
for the Asplund metallicity thus appear to pose similar difficulties for our understanding of 
M67 as for helioseismology, in the sense that stellar models with higher values of Z Q are 
better able to explain the observations. 

As the CNO elements play such a key role in the interpretation of the cluster data, 
one must first consider whether their abundances, relative to iron, are higher in M 67 than 
in the Sun. An increase in the abundances of C, N, and O by ~ 0.1 dex (or more, if the 
increase is limited to one or two of these elements) may be enough to reconcile the pre- 
dicted and observed CMDs. However, while a true cosmic spread among solar metallicity 



disk stars with an amplitude of a bout 0.1 dex in [C /Fel , [N/Fe], and 



excluded (see the recent surveys by lReddy et al.ll2003l and lEcuvillon et al 



O/Fej cannot be 



2006), the abun- 



dance analysis of M 67 main-sequence stars carried out by iRandich et al.l (120061 ) does not 
support the possibility that M 67 has [CNO/Fe] <1 0.1. These authors derived a mean oxygen 
abundance corresponding to [O/Fe] = —0.01 ± 0.03, which agrees well with the findings of 



1 Although one might question the reality of the gap in the Montgomery et al. (1993) CMD, there is no 
doubt about its e xistence (in the same magnitude range as that indicated in Fig. [2|) in the CMD published by 



Sandquistl (j2004T ) for the high-probability single-star members. Moreover, the observed distribution of single 



stars along the cluster fiducial sequence (see San dquist's Fig. 14) is ve ry reminiscent of that predicted by 
best- fitting isochrones (see Fig. 14 in the study by Michaud et al. 2004). [The CMD obtained by Sandquist 
is not used here because he did not provide BV photometry (only VI), an d his V — I colors, but not hi s 
V magnitudes, are at odds with other measurements (see the discussion bv I VandenBerg fc Stetson! 120041 ). 
Perhaps the main advantage of using BV data in our analysis is that, as mentioned above, the (B — V)-T e g 
relations that have been used to transpose the models to the observed plane appear to be especially reliable.] 
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Tautvaisiene et al.l (120001 ) for cluster giants. The latter also derived C and N abundances, 
though, in the case of giant stars, they are expected to have been altered by CN burning and 
the first dredge-up. They find that [C/H] and [N/H] are typically —0.2 and +0.2, respec- 
tively, which is roughly con s istent with initial C and N abundances slightly lower than solar. 
In fact, iFriel fe Boesgaardl dl992h found a m e an [C /H] of -0.09 ± 0.03 for three F dwarfs 
of M67. [Based on the work of iNissen et al.l (120021 ). we find that the expected differential 
effects due to the use of 3-D model atmospheres in the determination of [O /Fe] in the cluster 
turnoff stars, like those analyzed by Randich et al., would amount to +0.03 dex relative to 
the solar value.] Thus, it seems improbable that the CNO abundances relative to iron are 
enhanced in the cluster with respect to the Sun by as much as 0.1 dex. We also note that an 
increase of the Ne abundance, which has been suggested (by lBahcall. Basu fe Serenellill2005l ) 
to solve the discrepancy between solar oscillation data and the Asplund et al. abundances 
would not resolve the M 67 conundrum. 

Although our present understanding of diffusive processes precludes the possibility that 
there is a larg e differ ence between the interior and photospheric metallicities of the Sun, 



Michaud et al.l (12004J ) have found from their calculations for Z & > 0.0175 that the lowest 
mass diffusive model with a convective core on the MS is ~ 1.1OA4 , as opposed to 1.14A4© 
in the case of non-diffusive models. The difference in mass is comparable to the difference 
in A^tr predicted by our calculations when the assumed value of Z is increased from 0.0125 
(Asplund et al.) to 0.0165 (GS98). Consequently, it seems quite possible that \ow-Z diffusive 
models (with or without some core overshooting) may be able to satisfy the M67 constraint. 
Further investigation is needed. However, even if this leads to a consistent explanation for 
both the Asplund et al. solar abundances and the M67 CMD, the large discrepancies with 
helioseismology would remain. 



4. Conclusions 

The new Asplund et al. metallicity for the Sun presents some difficulties for fits of solar 
abundance models to the M 67 CMD, in that they do not predict a gap near the turnoff 
where one is observed. Whether or not stellar models predict a gap at the observed My 
in this open cluster depends quite sensitively on the assumed CNO abundances, and it is 
mainly the revision of these elemental abundances (along with that of Ne) that is responsible 
for the decrease in Z Q from 0.0165 (GS98) to 0.0125. Isochrones for the higher value of Z 
are able to reproduce the detailed CMD morphology of M 67 in the vicinity of the turnoff 
without apparent problems. Interestingly, it is primarily the reduction in the abundances of 
CNO that has resulted in substantial difficulties for helioseismology; e.g., an increase from 
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~ 0.3% to ~ 3% in the differences between the predicted and inferred sound speed squared 
profiles. Are solar oscillation studies and our investigation of M 67 telling us that the low 
solar Z determined by Asplund et al. is wrong? 



Not necessarily. It is possible, judging from the work of iMichaud et al.l (120041 ). that 
models which take diffusive processes into account may not have the same difficulties as 
the models used in this study, which neglect such processes. That is, if the Asplund et 
al. solar abundances are correct, only those \ow-Z models that treat diffusion may be able to 
reproduce the M 67 CMD. This possibility, which would not resolve the quandary presented 
by solar oscillations, needs to be carefully studied. Diffusion clearly adds another level of 
complexity to the problem since, e.g., the initial abundances of the gas out of which the Sun 
and M67 formed would have been somewhat different, if they have the same abundances 
today, because the Sun is up to 1 Gyr older than M67. Importantly, MS stars in M67 
should show systematic variations in their surface metallicities as a function of T e g due to 
the operation of diffusive processes (see Michaud et al.). 
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Fig. 1. — Grids of evolutionary tracks for the indicated values of Z and masses (in solar 
units). The sequences for 1.155 and 1.156A^ (dashed curves) and for 1.195 and 1.196j\4© 
(solid curves) show how rapidly the morphology of the tracks changes near _M tr (see the 
text). Although not plotted, both grids were extended down to 0.6M®. 
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Fig. 2. — Comparison of isochrones for the indicated and ages and metallicities with the 
CMD of M67 (Montgomery et al. 1993), on the assumption of the reddening and apparent 
distance modulus specified in the lower left-hand corner. The vertical line bounded by short 
horizontal lines, located adjacent to the turnoff, indicates our estimate of the luminosity 
spanned by the gap in the observed CMD. Only the dashed isochrone has a blueward hook 
in the same magnitude range. 
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Table I. Solar Elemental Abundances 



log A 



Element 


Grevesse & Sauval (1998) 


Asplund et al. a 


H 


12.00 


12.00 


Hc b 


10.9738 


10.9487 


C 


8.52 


8.41 


N 


7.92 


7.80 





8.83 


8.66 


Nc 


8.08 


7.84 


Ma 

in a 


O.oo 


D.oo 


Mg 


7.58 


7.58 


Al 


6.47 


6.47 


Si 


7.55 


7.51 


P 


5.45 


5.45 


s 


7.33 


7.33 


CI 


5.50 


5.50 


Ar 


6.40 


6.18 


K 


5.12 


5.12 


Ca 


6.36 


6.36 


Ti 


5.02 


5.02 


Cr 


5.67 


5.67 


Mn 


5.39 


5.39 


Fe 


7.50 


7.45 


Ni 


6.25 


6.25 


Y 


0.2676 


0.2559 


Z 


0.0165 


0.0125 



a Abundances for C, N, O, Ne, Si, Ar, and Fe were pro- 
vided by M. Asplund (2004, priv. comm.): Grevesse & 
Sauval (1998) abundances are assumed for all other ele- 
ments heavier than helium. 



b Dctcrmincd from Standard Solar Models. 



